Self-consistent supercell approach to alloys with local environment effects 



O. E. Peil 

/. Institut fiir Theoretische Physik, Hamburg University, Hamburg, Germany 

A. V. Ruban 

Royal Institute of Technology, Department of Material Science & Engineering, SE- 10044! Stockholm, Sweden 

B. Johansson 

Royal Institute of Technology, Department of Material Science & Engineering, SE-10044! Stockholm, Sweden and 
Condensed Matter Theory Group, Department of Physics and Materials Science, Uppsala University, SE-75121 Uppsala, Sweden 

We present an efficient and accurate mettiod for caicuiating eiectronic structure and reiated prop- 
erties of random alioys witti a proper treatment of local environment effects. Tlie method is a 
generalization of the locally self-consistent Green's function (LSGF) technique for the exact muffin- 
tin orbital (EMTO) method. An alloy system in the calculations is represented by a supercell with 
a certain set of atomic distribution correlation functions. The Green's function for each atom in the 
supercell is obtained by embedding the cluster of neighboring atoms lying within a local interac- 
tion zone (LIZ) into an effective medium and solving the cluster Dyson equation exactly. The key 
ingredients of the method are locality, which makes it linearly scaling with the number of atoms 
in the supercell, and coherent-potential self-consistency of the effective medium, which results in a 
fast convergence of the electronic structure with respect to the LIZ size. To test the performance 
and accuracy of the method, we apply it to two systems: Fe-rich bcc-FeCr random alloy with and 
without a short-range order, and a Cr-impurity on the Fe surface. Both cases demonstrate the 
importance of taking into account the local environment effects for correct description of magnetic 
and bulk properties. 



I. INTRODUCTION 

Accurate treatment of the electronic structure of dis- 
ordered systems is a highly non-trivial problem, which 
requires the use of a proper statistical model. In the 
case of metallic random alloys on an ideal crystalline 
lattice, the simplest statistical averaging can be done 
for a single site leading to the translationally invariant 
effective medium best given by the so-called coherent- 
potential approximation (CPA) .1^^ The CPA constitutes 
the basis of most of the first-principles techniques for the 
electronic-structure calculations of random alloys, where 
it is usually combined with the multiple scattering the- 
ory, or Korringa-Kohn-Rostocker (KKR) and related 
methods.EKS] Although the CPA has proven to be a quite 
successful approach for the electronic structure of many 
alloy systems, its application range is restricted either to 
completely random alloys or/and to alloys where local 
environment effects in the electronic structure are small, 
in the sense that their average is accurately given by the 
corresponding CPA effective medium. 

The proper account of the local environment effects 
in the electronic structure calculations, brought about 
either by correlated atomic distribution of alloy compo- 
nents (atomic short range order) or by multisite elec- 
tronic structure correlations, can be done in several 
ways which might be loosely subdivided into analyti- 
cal Green's-function-based approaches, in which the non- 
locality is taken into account by the summation of some 
subset of diagrams, and into supercell methods that make 
use of the self-averaging property of certain quantities, 
implying that the configurational averaging can be re- 



placed with the averaging over a sufficiently large super- 
cell. 

The main objective of analytical approaches is to ob- 
tain a full Green's function or, alternatively, a self-energy 
of a disordered system. This is achieved by summing 
a large subset (ideally all) of diagrams corresponding 
to multiple scattering from different sites either for the 
whole system, like it is done in the augmented space 
formalism,!^ or for a specific set of clusters represent- 
ing a given alloy system, as it is implemented in the 
cluster extensions of CPA.fi^HIIl The most developed ex- 
ample of this family of methods is the non-local CPA, 
that has been recently combined with first-principles 
calculations.'i^HIll this technique, the averaging is per- 
formed over all possible configurations of a specifically 
chosen cluster which tiles the entire underlying crystal 
lattice. In this case, as in similar analytical approaches, 
the average translational symmetry of the underlying lat- 
tice is lost together with the corresponding simple re- 
ciprocal space formalism, although it can be recovered 
with some additional effort. ^-I However, to our knowl- 
edge, there has been no implementation of the non-local 
CPA within density functional theory (DFT) strictly con- 
sistent with the total charge density. 

The requirement of the charge self-consistency is im- 
portant not only for accurate total-energy calculations, 
but also because the correct account of the local en- 
vironment effects themselves demands taking into con- 
sideration the response of the electronic density to the 
whole surrounding system, which can be fulfilled if the 
one-electron potential of every atom is determined in the 
DFT-self-consistent way. The difficulty of achieving this 
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self-consistency in Hamiltonian-based approaches, such 
as augmented space method, seems to be the main hur- 
dle on the way to establish them as an accurate tool for 
the first-principles calculations of disordered alloys. 

One of the ways to treat a disordered system within 
DFT self-consistently is to replace the configurational 
averaging by the averaging over a large sample, which 
is at the heart of supercell-based methods. These aver- 
agings are equivalent for an infinite-size sample and for 
quantities that possess the property of self-averaging, i.e. 
almost all observables of interest, including the Green's 
function and total energy. Moreover, if the interatomic 
interactions in the system are screened and hence short- 
ranged, as it is the case for most metallic alloys, the size 
of the sample can be chosen to be reasonably small with- 
out introducing a significant error.^^S 

In a naive, or direct, implementation, one would set up 
a supercell with periodic boundary conditions, solve it ex- 
actly by, for instance, evaluating the Green's function of 
the entire supercell within a DFT code, and then average 
the on-site components of the Green's function for every 
component to get the alloy Green's function. This ap- 
proach suffers from at least three serious problems: first, 
it can be computationally extremely demanding to find 
a Green's function of a large supercell, since the com- 
plexity usually scales as the third power of the supercell 
volume; second, the Bloch states of a periodic system, 
with their infinite lifetimes, cannot properly represent 
the states of a disordered alloy, and even after averaging 
over sites, the alloy Green's function would "feel" the pe- 
riodicity, which may result, for example, in an ill-defined 
conductivity; and finally, the Brillouin zone (BZ) of the 
underlying lattice is shrunk to a small BZ of the super- 
cell, and additional approximate routines are needed to 
restore the original translational symmetry. Note, that 
the problem with Bloch states can be partially resolved 
by performing additional averaging over different real- 
izations of the supercell, but this approach is obviously 
computationally expensive. 

A method that overcomes all of the named problems 
is the lo cally self-consistent Green's function (LSGF) 
technique! ^^ ! ^'' ! The method, as it will be described be- 
low, combines the best features of both the analyt- 
ical and supercell-based approaches and in addition, 
it is manifestly charge-self-consistent. In a histori- 
cal perspective, it can be considered as a generaliza- 
tion and extension of the multishell, or embedded- 
cluster method (ECM) ^ and the locally self-consistent 
multiple-scattering (LSMS) method. -'^^ 

The main idea of LSGF is that the local (on-site) 
Green's function for each atom in the supercell is de- 
termined from the Dyson equation restricted to the local 
interaction zone (LIZ) consisting of the given atom and 
its local environment embedded in the CPA-like effec- 
tive medium which is, in turn, built upon all the atoms 
residing on the alloy (sub)lattice. In this case, the peri- 
odic boundary conditions for the whole supercell are used 
only to determine the local environment of atoms close to 



the boundary of the supercell and for solving the Poisson 
equation in the electrostatic problem. The translational 
symmetry of the underlying lattice is used for solving the 
CPA equation for the effective medium and this method, 
thus, allows the use of the proper reciprocal-lattice for- 
malism, avoiding at the same time the ideal non-decaying 
Bloch states in the case of alloys, which plague the di- 
rect supercell approach. Besides the mentioned advan- 
tages, the computational complexity of LSGF scales lin- 
early with the supercell size. 

Certainly, the main playground of the LSGF method 
is the consistent and accurate calculations of the elec- 
tronic structure of random alloys with or without short 
range order. Also, increasing the size of the LIZ allows 
the systematic investigation of the local environment ef- 
fects in alloy systems. Besides, it is a powerful tool 
for ev aluat ing screened Coulomb interactions in random 
alloys .'^Sl^i] At the same time, the application range of 
the method reaches far beyond homogeneous disordered 
systems. For example, by judiciously choosing the size 
of the LIZ to hide periodic images in a slab geometry, 
one can efficiently treat problems involving impurities on 
or near surfaces. In a similar way, the explicit evalu- 
ation of the interatomic interactions between a pair or 
multiple of impurities in the host metal is possible. The 
Green's-function approach embodied in the method pro- 
vides additional benefits in treating paramagnetic ran- 
dom systems. Some examples of such applications will 
be presented in this paper. 

The previous practical realization of the LSGF method 
was done using the KKR m ethod within atomic sphere 
approximation (ASA).Elll The KKR- AS A method itself 
suffers from the normalization error: electronic states 
are normalized within Wigner-Seitz spheres instead of 
Wigner-Seitz cells. This leads, in turn, to a quite sub- 
stantial error in the total energy, making, for instance, 
practically impossible the accurate calculations of the 
total-energy variations related to the deformation of the 
crystal structure. At the same time, the LSGF approach 
is a general formalism that can be implemented in any 
code based on multiple-scattering theory. In this paper, 
we describe the implementation within the exact-muffin- 
tin-orbital (EMTO) method.i^Stsn! The method belongs 
to the family of screened-KKR techniques with the ba- 
sis set formed by so-called third-generation muffin-tin 
(MT) orbitals introduced by Andersen.I^S During the last 
decade, this method, combined with the fuU-charged- 
density technique, has proven to be sufficiently accurate 
in calculations of a wide spectrum of alloy properties.!^ 

The paper consists of two main parts. In the first 
part, namely Sectionjllj we describe the LSGF formalism, 
starting with some details of the EMTO method that are 
essential for understanding the implementation. In the 
second part. Sect ion [Till we provide the results of test cal- 
culations followed by real calculations demonstrating the 
ability of the method to treat both homogeneous random 
alloys, magnetic, as well as non-magnetic, and inhomo- 
geneous systems, such as surfaces. In the last section, we 
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conclude the results and briefly discuss possible exten- 
sions of tlie current LSGF implementation. 



II. METHOD 

As has been mentioned in the Introduction, in LSGF, 
one starts from building a supercell with periodic bound- 
ary conditions by replicating a unit cell of the underly- 
ing lattice. The supercell is populated with components 
in any desired (random or ordered) configuration (see 
IIB| . The main ingredients of the EMTO- 



Subsection 

LSGF (ELSGF) approach are the supercell setup, EMTO 
method itself, and LSGF part involving the solution of 
the restricted Dyson equation. 

The EMTO part (more generally, the KKR part) of 
the ELSGF method runs similar to the usual EMTO- 
GPA implementation, as far as it concerns the effective 
medium, density of states, charge density, and total en- 
ergy. The full symmetry of the underlying lattice is em- 
ployed, reflecting the single-site character of the effective 
medium. After the effective medium is found, the cluster 
path operator is evaluated using the ECM, or restricted 
Dyson equation, for each atom in the supercell. This 
cluster path operator is subsequently used for the nor- 
malization of states, determining the density of states, 
charge density, and total energy. 

In this section, wc first briefly outline the EMTO 
Green's-function formalism and the features specific to 
its implementation within the LSGF technique. Then, 
we present strategies for choosing a proper supercell, and 
finally describe the LSGF formalism within EMTO. De- 
tails concerning the structure constants, the construction 
of the optimized one-electron potential, the full charge 
density, as well as the total energy calculations within 
the full-charge-density formalism, are the same as in the 
usual EMTO method, and their comprehensive descrip- 
tion can be found in Ref. [23 and EHl 



A. The EMTO Green's-function formalism 

The main idea behind the EMTO method is to keep 
the simplicity of the MT- or screened-KKR method by 
making use of the muffin-tin geometry of the one-electron 
potential, but to improve the accuracy of the electronic 
structure calculations to the level of the full potential 
methods. The latter can be achieved by replacing usual 
MT-spheres with large overlapping MT-spheres which al- 
low a better approximation for the full potential. At the 
same time, additional non-overlapping screening spheres 
are used to define boundary conditions for the solutions 
(referred to as screened spherical waves) in the interstitial 
region. 

The path operator, g{z), is defined as follows. 
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with the kink operator 

KwL'RLiz) - aR'S'R,L'RL('«^) - Srr' Sll' aRD^i{z) , 

(2) 

where or, is the screening-sphere radius at site R, 
^r'L'rlI'^^) screened structure constants depending 
on energy = z — vq defined with respect to the muffin- 
tin zero vq, and D^i{z) the potential function determined 
as usual by the logarithmic derivative of the partial waves 
at the MT sphere. Note the sign convention for the po- 
tential function and structure constants. 

For systems with translational symmetry, the on-site 
path operator is determined as 
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where integration is performed over the Brillouin zone 
of the crystal lattice. We write all expressions for a one- 
atom Bravais lattice for clarity; the extension to multiple- 
atom basis is straightforward. 

To evaluate the density of states (DOS), the path op- 
erator must be properly normalized with the overlap ma- 
trix. The overlap matrix in the EMTO formalism is given 
by ^r.'l'rl(-2)j the energy derivative of the kink opera- 
tor. The number of states is then, 

^f{s) - - ^ Jdz [G{z) - Dp^U^)] , (4) 
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where the integration along the half of the contour em- 
bracing the valence band below energy e is performed, 
and pole contributions, -Dpoies: coming from the poles of 
1/D{z) and D{z), are subtracted. The Fermi energy, ep, 
is found from the condition Af{ep) = N^i, where N^i is 
the number of the valence electrons. 

In the case of a random alloy on a lattice, the CPA 
equations are used to determine the electronic structure 
given by the coherent path operator, g, through the cor- 
responding coherent potential operator, D, of the single- 
site effective medium: 
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The path operators of the i-th alloy component, gi, are 
found via the single-site Dyson equation, 

9i ^ 9 + 9{D - Di)g,, (7) 

from which the coherent path operator is determined as 

E " 9, (8) 
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where Ci is the concentration of alloy components. The 
last three CPA non-linear equations are solved self- 
consistently. 

In EMTO-CPA, the correctly normalized Green's func- 
tion and number of states (per Wigner-Seitz cell) are de- 
termined as 
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and the pole contributions are weighted by concentra- 
tions Ci accordingly. 



B. Supercell 

In the LSGF calculations, an alloy system is repre- 
sented by a supercell model. In general, creating a su- 
percell with needed statistical properties, given by its 
atomic-distribution correlation functions, is a highly non- 
trivial task, mathematically equivalent to the optimiza- 
tion of a many-variable function in a multidimensional 
space. The initial building block of the supercell is de- 
termined by the underlying lattice containing Nq basis 
atoms. Let us note that the choice of the initial unit cell 
of the supercell is quite arbitrary, and such a unit cell can 
be different from that of the underlying crystal lattice in 
the subsequent LSGF calculations, provided that the su- 
percell is conformal to the underlying lattice. A simple 
example is the choice of the cubic unit cell containing two 
atoms in the case of bcc structure as a building block of 
the supercell and the use of the bcc translational sym- 
metry in the LSGF calculations. In general, the choice 
of the unit cell is motivated by the model of an alloy 
system that may have several different sublattices with 
different compositions and distributions of alloying ele- 
ments (for an application of ELSGF to a rather complex 
example of the FeCr cr-phase with 30 atoms per unit cell 
subdivided into 5 nonequivalent sublattices with distinct 
compositions, seeP3). 

The setup of the supercell starts from the definition of 
a desired alloy configuration by characterizing its atomic 
correlation functions that can be defined in different 
ways. For a homogeneous binary alloy, for instance, they 
can be given by the average products of spin-like vari- 
ables, (Ti, taking on values -1-1 or —1, depending on which 
alloy component occupies site i: 
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where ^y"'* is the n-site correlation function for cluster /, 
and (...) is the average over the supercell. A completely 
disordered configuration is given by ^jr"-* = cr", with a = 

2c — 1 = ^^^^ where c is the concentration of one of the 
alloy components. 



Although this definition is easily generalized to the case 
of multicomponent and inhomogeneous random alloys, 
where different sublattices with different alloy composi- 
tions and atomic short range order are present, it is ex- 
tremely difficult to use it in practice. The main obstacle 
here is the finite and quite restricted size of the super- 
cell. For instance, the number of sites in the supercell, 
N, defines the possible concentrations to be only k/N, 
where k = 0,1, ... N. Much more severe restrictions to 
the possible pair and multi-correlation functions origi- 
nate from the geometry of the underlying lattice, e.g., 
from the number of possible clusters (coordination num- 
ber for pair correlation functions). Nevertheless, a su- 
percell consisting of about 500-2000 sites is usually suf- 
ficiently big to model a large variety of alloy systems, 
especially taking into consideration our restricted knowl- 
edge about atomic distribution correlation functions in 
real alloys. 

The main condition for a supercell to be a valid rep- 
resentation of a random alloy with a specific short range 
order is to have the same correlation functions as the 
given alloy for those clusters, /, (or coordination shells 
in case of pair correlation functions) which affect a phys- 
ical property of interest. If the cluster expansion is ap- 
plicable to a particular observable, A, then the relevant 
clusters are those, for which expansion coefficients, Af, 
have non-zero (in practice non-negligible) values: 
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Here, we do not discuss for which physical properties 
such an expansion can be valid in general, and how fast 
the expansion converges, but as an example we can men- 
tion the total energy of an alloy which, as usually tacitly 
assumed, can be expanded in this way. The coefficients 
Af are then just the effective interactions of the corre- 
sponding Ising Hamiltonian. 

Let us note that a finite supercell of a restricted size 
cannot represent a random alloy in general (and called 
"random" for that matter): the same supercell can be 
"random" for the same system for one property and "or- 
dered" for the other, not mentioning different systems. 
In principle, one should check the configurational depen- 
dence of the observable of interest (or ideally to find out 
coefficients Af) prior to using a supercell as a random- 
alloy model. For instance, the expansion coefficients of 
the total energy are related to the effective interatomic 
interactions that can be found either with the aid of the 
generalized perturbation method within single-site CPA 
or with the cluster inversion method. 

In LSGF, the real correlations of the supercell are only 
taken into consideration inside the LIZ, while outside the 
LIZ, the correlation functions correspond to a completely 
random alloy. An observable calculated by the LSGF 
method is given by the following formula: 
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cr" are the correlation func- 



where Cf = Q 
tions of the completely random alloy. The condition for 
the cluster to belong to the LIZ is that one of its vertices 
coincides with the position of the central atom, and all 
cluster atoms belong to the LIZ. This condition will be 
elaborated in the next sections, where we give the details 
of the electronic structure calculations with the LSGF 
method. 

The electrostatic energy is calculated "exactly" for a 
given supercell. This means that the summation is not 
restricted to the LIZ but is rather performed over the 
entire periodic infinite system. In principle, this may 
result in a spurious electrostatic interaction between pe- 
riodic images. Fortunately, in most metallic alloys, the 
pair interactions are screened, the screening length being 
rather small. The important condition that the range 
of the pair-correlation function and the screening length 
are within a supercell can, thus, be easily satisfied for 
supercells containing several hundred atoms. 
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multiple-scattering problem for the LIZ-cluster embed- 
ded into the effective medium, which enables one to take 
the local environment effects into consideration. The size 
of the LIZ is usually defined as a number of coordination 
shells constituting the cluster centered around a given 
atom, with the LIZ of size one (LIZ=1) corresponding to 
the cluster consisting of a single central atom. The rest 
of the system outside the LIZ is given by the effective 
medium. 

The multiple-scattering problem for such a setup is 
solved exactly with the aid of the Dyson equation as it is 
done in the Embedded Cluster Method (ECM).ISl Con- 
sider a LIZ-cluster given by atoms with potential func- 
tions Di, where i runs over cluster sites R^. The cluster 
is embedded into the effective medium defined by a real- 
space path operator gij and a coherent potential function 
Di. The path operator of the cluster, gij, is then found 
within the ECM as 



9ij — 9ij 
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(18) 



from which one finds immediately 



Given a supercell containing N atoms, the calculation 
of the electronic structure within the ELSGF method 
starts from the determination of the translationally in- 
variant CPA effective medium built upon all the atoms 
on a corresponding (sub)lattice. In the simplest case of a 
Bravais lattice, corresponding CPA equations arc similar 
to Eqs. (|6f-(|8|: 
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(15) 
(16) 



Here, gQ and D are the on-site effective-medium path op- 
erator and logarithmic derivative; gi is the on-site path 
operator of site i. These nonlinear equations are solved 
self-consistently for D and g for a given set of one-electron 
potentials in the supercell. Eqs. (14|-(16) are the CPA 



equations for an iV-component alloy. This makes the ef- 
fective medium single-site-self-consistent and all transla- 
tional and point symmetries of the underlying lattice are 
preserved. In addition, a connection to the CPA guaran- 
tees the analyticity of the effective medium and hence of 
its real-space path operator determined as 
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where R,;, Rj are the positions of sites i and j of the 
lattice. 

Once the effective medium is defined, the electronic 
structure for every site can be determined by solving the 
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The Fermi energy of the system, ep, is determined 
from the normalization condition, A/'(£f) = A^e^ for the 
number of electron states in the supercell, N^i, where the 
number of states is defined in a way similar to that in the 
EMTO method: 



AA(£) = -^ J dz [G(z) - i?poie,(z)] 



(20) 



with the Green's function G{z) having an additional con- 
tribution compared to that of the EMTO-CPA Green's 



function (Eq. (lOl). The point is that unlike the EMTO- 



CPA, where the states are normalized within a single-site 
Wigner-Seitz cell, the states in ELSGF are normalized 
within the entire LIZ. Starting from the expression for 
the Green's function in EMTO-CPA (10), rewriting it in 



the real space, and replacing the effective-medium path 
operator by the cluster path operator within the LIZ, one 
gets 
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dk 
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The electron density for each site is obtained from the 
on-site path operator, gi =gii, in the same way as it is 
done in the EMTO method.I^^ The DFT self-consistency 
loop is then closed by evaluating the one-electron poten- 
tial for every site, assuming the translational symmetry 
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(periodic boundary conditions) for the entire supercell in 
order to solve the electrostatic problem exactly for each 
site. Finally, after the self-consistency is reached, the full 
charge density is determined in order to perform accu rate 
calculations of the total energy of the supercellP^^H 



D. Disordered-local-moment model for the LIZ 

Accurate calculations of the electronic structure and 
energetics of paramagnetic alloys with the finite magni- 
tudes of local magnetic moments on atoms present a chal- 
lenge for modern first-principles methods. As has been 
proven in Ref. if the magnitude of the local magnetic 
moments does not fluctuate strongly and the spin-orbit 
coupling is negligible, such a state is accurately described 
by the collinear disordered-local-moment (DLM) model, 
where atoms with spin-up and spin-down orientations of 
their local magnetic moment are distributed randomly 
on the underlying lattice. 

Although the DLM configuration can be modeled by a 
supercell with randomly distributed atoms having differ- 
ent spin orientations, such a naive supercell representa- 
tion can, in fact, lead to incorrect results because in re- 
ality, magnetic degrees of freedom fluctuate (transverse 
fluctuations are implied here) rapidly and create, thus, 
a local environment different from the one with static 
magnetic moments. In this sense, a CPA-based scheme 
seems to be a better method to calculate the systems 
in the DLM state. Besides, a supercell model with ran- 
domly distributed static moments becomes too cumber- 
some, since specific atomic correlations should be set up 
not only between alloy components, but also between 
their spin-up and spin-down counterparts. 

Within the LSGF method, a straightforward DLM- 
CPA implementation is possible only in the single-site 
mode (LIZ=1), when the required CPA averaging and 
thus the Dyson equation for every site and spin can 
be solved for the appropriately spin-averaged effective 
medium. However, this simple scheme obviously breaks 
down when the nearest-neighbor atoms are included in 
the LIZ. The problem here is that the potential functions 
of the neighbors of the central atom would incorrectly 
correspond to a specific magnetic configuration rather 
than the random one, as it should be in the DLM state. 

Clearly, a correct description of the DLM state implies 
here that the central atom " sees" its neighbors inside the 
LIZ in spin-averaged states. Such a solution can be ef- 
ficiently implemented in the LSGF method by perform- 
ing partial constrained averaging of the spin states for 
all the sites inside the LIZ except the central one. This 
amounts to choosing an appropriate effective potential of 
a given atom in the same way as it is usually done within 
CPA-DLM. To be more specific, we start by defining the 
on-site path operators in two spin-channels by solving 



corresponding single-site Dyson equations, 
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A partially averaged (DLM-averaged) path operator is 
then introduced for each site and global spin channel a 
as 
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The paramagnetic effective medium is now determined 
from these DLM-averaged path operators, with the self- 
consistency condition being 
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Once the effective medium is determined, one can use 
the following equation 
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to find a corresponding DLM-averaged potential function 
for each site, 



(27) 



We then solve Eq. ( 18 ) for each site assuming that all 



the atoms in the LIZ surrounding the central one are 
in the paramagnetic state. The potential function for 
each non-central atom in the cluster is thus replaced by 
a DLM-averaged potential function. Specifically, for a 
given cluster with a central site i, the cluster potential 



function Z?, is defined as follows, 



3l G 



3=1': 
otherwise. 



(28) 



III. TESTS AND RESULTS 

A basic test for the ELSGF method is the conver- 
gence of quantities of interest (the total energy in the 
first place) with respect to the size of the LIZ. By con- 
struction, the method observes two limits: LIZ — cxi 
[N — > oo), corresponding to the formally exact solution 
of the Dyson equation for the entire system; LIZ = 1, 
equivalent to the CPA with the correct account of electro- 
statics (sometimes referred to as the isomorphous CPA). 
In between these two limiting cases, the convergence with 
respect to the LIZ depends pretty much on the observ- 
able or, more strictly, on how fast the cluster expansion 
coefficients decay with distance, as has already been em- 
phasized in Sec. |IIB[ In particular, the expansion coeffi- 
cients of the total energy are related to effective interac- 
tions, and the convergence test can be used as a rough 



estimate of the range of the effective interactions!^ An 
example of such a calculation is given below. 

Later in this section, we demonstrate some of the ca- 
pabilities of the ELSGF method by applying it to real 
systems. Emphasis is made on the effects of short-range 
order, especially in magnetic systems. Also, a surface- 
segregation problem is considered as an example of an 
inhomogeneous system. 



The range of the relevant correlation functions 
from ELSGF 



The direct relation between the convergence with the 
LIZ size and the range of the correlation functions makes 
it possible to estimate the latter by varying the size of 
the LIZ in LSGF calculations. As has been discussed 
in Sec. |IIB[ if an observable is self-averaging and can 
be expanded in terms of the atomic-distribution correla- 
tion functions, the expansion is given by Eq. ( 13 1 within 



LSGF. According to this equation, the correlation func- 
tions beyond the LIZ correspond to those of a completely 
disordered alloy, and thus, by calculating a completely or- 
dered alloy with LSGF, one captures only contributions 
from the correlation functions corresponding to the LIZ- 
cluster. In view of this, one can estimate the range of the 
relevant correlations as the minimal size of the LIZ that 
provides the same result as the one given by a direct ab 
initio calculation for the given ordered structure. 

As an example, we present here results for the DOS 
and the total energy of a completely ordered B2-NiAl 
phase, also considered in Ref. [T71 Let us note that this 
ordered phase represents the worst possible case for the 
LSGF method since its every coordination shell consists 
of atoms of only one type while in a random alloy with 
the same equiatomic composition, the number of atoms of 
both types should be equal (on average). Formally, this 
is given by the corresponding atomic correlation func- 
tions, or Warren-Cowley SRO parameters, which take on 
values ai = —1, 1, 1, —1, 1, 1, —1, . . . for the first several 
coordination shells in the B2 structure. Here, —1(1) cor- 
responds to the case when every atom has only atoms of 
the opposite (same) type at the corresponding coordina- 
tion shell, while the SRO parameters are zero in the case 
of a random alloy without short-range order effects. 

In Fig. [ij we show the DOS of the B2-NiAl phase ob- 
tained in the ELSGF calculations with different sizes of 
the LIZ, together with the DOS calculated by the EMTO 
method. One can see that at least four coordination 
shells (LIZ=5) need to be included in the LIZ to repro- 
duce the main features of the DOS with a reasonable 
accuracy, and still some small DOS features are not well 
reproduced even when eight coordination shells (LIZ=9) 
are included. 

A similar but a bit more interesting example is shown 
in Fig. [2] where the ELSGF calculations of the DOS of 
B2-FeSi presented. This system has unusual magnetic 
properties very sensitive to the local environment. In 
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FIG. 1. (Color online) The DOS for B2-NiAI obtained in 
the ELSGF calculations with different LIZ sizes (thick red) 
compared to a direct calculation (shaded region). LIZ = 1 
corresponds to a single-site approximation but different from 
CPA, because electrostatic potential is calculated for the or- 
dered structure. 



particular, there is quite large magnetic moment in bcc 
random Fco.sSio.s alloy, but it disappears in the B2 or- 
dered structure. Performing the ELSGF calculations of 
the B2 structure with increasing LIZ one can find the ef- 
fect of every coordination shell on the magnetic moment. 
As one can see in Fig. [2) there is an obvious splitting 
of the bands in the single-site approximation (LIZ=1, 
equivalent to the isomorphous GPA, but with incorrect 
electrostatic potential, which is determined for the B2 
structure) . 

The inclusion of the first coordination shell in the LIZ 
(LIZ=2), results in a substantial reduction of the split- 
ting. In this case Fe atoms can "see" nearest neighbor 
Al atoms (and vice versa, Al atoms are surrounded be 
Fe atoms at the first coordination shell) while the rest 
of crystal is represented by the CPA effective medium. 
Thus, such a large reduction of the magnetic moment of 
Fe seems to be obvious in this case. The next, second, 
coordination shell of Fe atoms in the B2 structure con- 
sists of only Fe atoms, and thus one could expect a slight 
increase of the magnetic moment due to additional Fe- 
Fe interactions when the LIZ size increases from two to 
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FIG. 2. (Color online) The DOS for B2-FeSi obtained in the 
ELSGF calculations with different LIZ sizes (thick solid-red 
and dashed-cyan) compared to a direct calculation (shaded 
region). LIZ — 1 corresponds to a single-site approximation 
but different from CPA, because electrostatic potential is cal- 
culated for the ordered structure. 



three. However, the sphtting and magnetic moment are 
further reduced, and when the next, third coordination 
shell is included in calculations, the magnetic moment 
becomes practically zero (see Fig. [3|. This indicates, 
that most probably the Fe-Fe exchange interaction pa- 
rameters are antiferromagnetic in this case. 

For comparison we also show in Fig. |4] the DOS of ran- 
dom Feo.sSio.s alloy calculated by the ELSGF. In this 
case the supercell consisting of 256 atoms (4 x 4 x 8(x2)) 
atoms has been used, whose atomic distribution pair cor- 
relation functions were as in the random alloy up to the 
8th coordination shell. The DOS for LlZ—1 (isomor- 
phous CPA model) is actually very close to that for the 
B2 phase shown in Fig. [2] The inclusion of the first co- 
ordination shell in the LIZ leads to a slight modification 
of mostly the spin-majority band, while the inclusion of 
more distant coordination shells to a slight change of the 
minority spin-band. In Fig.[4]wc show only the result for 
LIZ=6 (five coordination shells included in the calcula- 
tions). However, this result is practically indistinguish- 
able from those for the LIZ=4 and 5. 

The relatively fast convergence of the DOS with the 
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FIG. 3. (Color online) Local (average over the random al- 
loy) magnetic moment of Fe in the ELSGF calculations of B2 
structure and an equiatomic random alloy as a function of the 
LIZ size. 



LIZ size for random alloys is a natural feature of the 
LSGF calculations with the CPA effective medium. The 
better the CPA works, the faster convergence. In most 
cases, the inclusion just of the first coordination shell in 
the LIZ provides very accurate description of the DOS 
of random alloys. In the case of Fe-Si it is a bit slow, 
exhibiting distant "local environment effect", which are 
most probably connected with the non-trivial magnetism 
in this system, also showing up in the ELSGF calcula- 
tions of the B2-FiSi. 

The evolution of the density of states with increasing 
LIZ gives an idea of how fast the electronic structure ap- 
proaches the one for an ordered system with the inclusion 
of the corresponding interactions, thereby providing an 
estimate of the range of the correlation functions respon- 
sible for specific features of the DOS. An important and 
interesting point here, however, is that although the DOS 
is directly related to the total energy (specifically, to the 
band energy), its details are not important for the en- 
ergy value, and the convergence of the total energy with 
respect to the LIZ size can thus be much faster. This is 
partly due to the integral dependence of the total energy 
on the DOS. This point is illustrated in Figs. [5] and |6j 
where we show the convergence of the total energy of B2 
NiAl and FeSi, respectively, as a function of the LIZ size. 

One can clearly see that in the case of the B2-NiAl, 
already starting from the LIZ corresponding to the first 
coordination shell (LIZ=2), the difference in the total 
energy becomes very small and remains so for larger LIZ. 
As was shown in Ref. \17\ the change of the total energy 
with the LIZ size can be traced back to the strength 
of the effective intersite interactions entering the Ising 
configurational Hamiltonian. For this particular case of 
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FIG. 4. (Color online) The DOS of the random bcc Feo.sSio.s 
alloy obtained in the ELSGF calculations with different LIZ 
sizes: LlZ—1 corresponds to the single-site CPA (shaded 
area), LIZ=2 (thick dashed red line) and 6 (thick solid blue 
line) are the results with local environment effects up to the 
first and fifth bcc coordination shells, respectively. Spin- 
majority band is shown in the upper panel and minority in 
the bottom panel. 



NiAl binary alloy, it can be written as 



(29) 



where summation runs over coordination shells, p, and 

the corresponding sites of the lattice, i and j; Vp are 
effective pair interactions and Sci = Ci ~ c the concen- 
tration fluctuation of the occupation number, q, which 
takes on values 1 and 0, depending on which alloy com- 
ponent occupies site i. 

For the system considered, one can, in fact, estimate 
the value of effective interactions. Fig. [5] shows that the 
strongest interaction is the one related to the first coordi- 
nation shell. Assuming that it is the pair effective inter- 
action that is dominant at the first coordination shell, its 
value could be assessed as just the difference between the 
total energy of a random alloy and that of the ordered one 
with LIZ=2 (which includes the first coordination shell). 



i.e. v}^' '-^ 32 mRy (factor 1/8 in Eq. (29) is canceled ex- 




FIG. 5. (Color online) Convergence of the total energy with 
respect to the size of the LIZ for B2-NiAl. The difference 
between ELSGF (EMTO-CPA - dotted hne) and EMTO total 
energies are shown. Red solid line (open circles) - MT-basis 
with /max ~ 3; blue dashed line (filled circles) - MT-basis 

with /max = 2. 
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(Color online) Convergence of the total energy with 
to the size of the LIZ for the B2 and random FeSi 



actly by the bcc coordination number 8). However, one 



should be aware of the fact that such a simple estima- 
tion is valid only in the case when the contribution from 
multisite interactions is small compared with that from 
the pair interactions. At the same time, there are quite a 
few multisite interactions contributing even in the case of 
LIZ=2, when only the first coordination shell is included 
in the calculations since they are for the all possible clus- 
ters in the LIZ with one site being the central site of the 
LIZ. 
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One can also see in Fig. [5] that the total energy of 
the B2-NiAl for LIZ=1, which is the single-site approx- 
imation, is not equal to that of the random alloy ob- 
tained from the CPA, ahhough LSGF with LIZ=1 is for- 
mally equivalent to the CPA. The reason for the differ- 
ence is the electrostatic energy, which is calculated in 
LSGF for the entire supercell representing an ordered 
system in our case. The difference between the LSGF 
result with LIZ=1 and the energy of the random alloy 
can be roughly associated with the electrostatic part of 
the effective interactions, or a so-called intersite screened 
Coulomb interaction.l^Sl 

In Fig. |6] we show the results of similar caculations 
but for the ordered-B2 and random FeSi alloys. One can 
see that the convergence of the total energy of the B2 
phase with respect to the LIZ size is worse than in the 
case of B2-NiAl. Another reason for that, apart from the 
convergence of the effective interactions, can be the fact 
that FM state is stable up to LIZ=3. The existence of 
quite long-range local environment effects in random FeSi 
alloy have been already discussed above for the DOS. It 
also leads to the unusually slow convergence of the total 
energy of random alloy, which is as a rule quite accurate 
already for the LIZ=2, when the first coordination shell 
is included in the LIZ. 



B. Short-range order effects in random FeCr alloys 

One of the advantages of the LSGF method is the fact 
that the computational complexity scales linearly with 
the number of atoms, which renders possible accurate 
first-principles calculations for systems containing up to 
several thousand atoms. Such a size of a supercell allows 
one to model a large variety of alloys with various con- 
centrations and atomic short-range orders (SRO). In this 
section, we demonstrate such a possibility for the case of 
ferromagnetic Fe-rich Fe-Cr alloys which have attracted 
great attention of scientists in different fields owing to 
potential applications of these alloys in industry. 

From a scientific point of view, this system is quite in- 
teresting and complicated, when it comes to its accurate 
first-principles description. One of the reasons is a com- 
plex interplay of magnetism at zero as well as at elevated 
temperatures, and interatomic interactions, and conse- 
quently, thermodynamic properties of these alloys .I^^HSll 
In particular, the type of alloying abruptl y changes with 
Cr concentration and temperature^^I^SEni in the compo- 
sition range of up to about 20 at.%Cr. While there is a 
quite strong ordering tendency between Fe and Cr atoms 
at low temperatures in the ferromagnetic state and at low 
Cr concentration, Fe-Cr alloys exhibit a phase separation 
behavior at higher temperatures, close to and above the 
Curie temperature, and with increasing Cr content. It 
was also demonstrated in Ref. 134| that the usual Ising 
model breaks down for this system due to a strong local 
environment dependence of the effective interactions. 

It is clear that real Fe-Cr alloys must have a certain 



amount of atomic short-range order, and its type and 
magnitude depends on the thermal treatment of alloy 
samples. At the same time, practically all calculations 
for random Fe-Cr alloys are done for completely random 
alloy configurations (in fact, with rare exceptions, these 
are just CPA-based calculations). In this section, we in- 
vestigate the effect of the atomic short-range order on 
some of the ground-state and elastic properties of the 
Fe-rich FeCr alloys using the ELSGF method. 

In Fig. [7j we show the generalized-gradient- 
approximation (GGAJP^ results obtained by the usual 
EMTO-CPA method and by the ELSGF for the depen- 
dence of the equilibrium theoretical lattice constant on 
Cr concentration and the (Warren-Cowley) SRO param- 
eter at the first coordination shell. 



The lattice constant is shown in relative units of the 
deviation from the average values (given by Vegard's 
law) for a given concentration.Esl The EMTO-CPA re- 
sults are similar to those obtained in Ref. [M] and the 
ELSGF results have been obtained using a 256-atom su- 
percell (8x4x4(x2) based on the cubic unit cell of the 
bcc lattice) for two alloy compositions of Feo.9375Cro.0625 
and Feo.875Cro.125. In the first case, calculations have 
been done only for a completely random alloy without the 
SRO, while in the latter case, we have calculated the equi- 
librium lattice spacing of four supercells with ai — —0.1, 
0, 0.1, and 0.2. All these supercells are, of course, not 
ideally random, but the deviation of the other correlation 
functions from those in a random alloy has been smalhl^Sl 

As one can see, the ELSGF results for the random alloy 
are below the corresponding CPA results, hence closer to 
the experimental values.^"* One of the reasons is that the 
EMTO-CPA self-consistent calculations have been per- 
formed for fixed values of screening constants entering the 
definition of the on-site screene d Cou lomb interactions in 
the single-site DFT formalism. I^^l2il Although they were 
determined for every composition,'^ they were kept con- 
stant in the total energy calculations for different lattice 
constants. Local environment effects can also play an im- 
portant role in this alloy, since they affect the magnitude 
of the magnetic moment of Cr atoms dramaticallyl^^^^^ll 
A clear manifestation of this effect is the dependence of 
the equilibrium lattice constant in Feo.875Cro.125 on the 
SRO parameter: it decreases proportionally to the value 
of ai (see also the top panel of Fig. |8]) 

Let us note that the result is counterintuitive, for neg- 
ative values of the SRO parameter mean ordering ten- 
dency, i.e. preferential occupation of the corresponding 
coordination shell by the atoms of the opposite type, 
while positive values correspond to phase separation. 
The reason for such a reversal is the behavior of the local 
magnetic moments of Cr atoms. Their magnitude in Fe- 
rich Fe-Cr alloys is roughly proportional to the number 
of Fe nearest neighbors, and they, thus, go down with 
increasing the SRO parameter at the first coordination 
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FIG. 7. (Color online) The dependence of the lattice constant 
on Cr concentration, ccr, and short-range order parameter on 
the first coordination shell, qi, from the EMTO-CPA (circles) 
and ELSGF (stars) calculations (the solid line is a guide to 
the eye). Deviation of the lattice constant from the Vegard's 
law is plotted. The only value obtained with ELSGF for ccr = 
6.25% is evaluated for a completely disordered alloy (a; — 0.0 
up to the 8th coordination shell). 



shell, as one can see at the bottom panel of Fig. [S] De- 
creasing magnetic moments of Cr results in the reduction 
of magnetic pressure, producing the effect of contraction. 
This example shows that SRO effects are especially im- 
portant in magnetic alloys. 



C. Elastic constants of Feo.875Cro.125 

The full-charge-density formalism implemented in the 
EMTO method allows one to perform relatively accurate 
calculations of elastic properties of solids.l^^ Although 
theoretical results are usually in a quite good agreement 
with available experimental data, the accuracy of such 
calculations, in particular related to the use of the singe- 
site CPA in the DFT self-consistency, is not known. 

When no alloys on sublattices are present, the ELSGF 
and EMTO methods are equivalent. In other words, the 
calculations of ordered systems are equivalent in accu- 
racy, and this makes possible the direct comparison of 
the results on elastic properties obtained by the EMTO- 
CPA and ELSGF methods. To this end, we have again 
chosen Fe-Cr alloys, whose elastic constants have been 
recently calculated by the EMTO-CPA method. In par- 
ticular, we have calculated shear elastic moduli, c' and 
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FIG. 8. (Color online) The dependence of the lattice constant 
(circles), bulk modulus (boxes), and magnetic moment of Cr 
(dots with whiskers) on the short-range order parameter on 
the first coordination shell, ai. The whiskers in the bottom 
panel show the mean-square-root of the distribution of the 
magnitude of magnetic moments on Cr atoms. 



C44, of the Feo.875Cro.125 alloy at the experimental lattice 
constant, 2.869 A, in both the ferromagnetic and param- 
agnetic, i.e. DLM, states. 

The ELSGF calculations have been done for a com- 
pletely random alloy and with some amount of atomic 
SRO at the first coordination shell, ai =0.1. The alloys 
have been modeled by 512-atom supercells (8x8x8 based 
on the primitive unit cell of the bcc lattice) . The usual 
EMTO-CPA calculations have been done exactly for the 
same set up of k-points (39x39x39 division of the full 
Brillouin zone in the Monkhorst-Pack methocP^l)^ Imax 
cut-off (=3 for the partial waves inside atomic spheres), 
and other parameters. 

For random alloys, an important difference between 
the EMTO-CPA and ELSGF calculations is that the 
DFT-based calculations of random alloys within single- 
site CPA should take into consideration the shift of 
the one-electron potential due to the on-site screened 
Coulomb interactions, V^^^r 120121] ^^^^ been demon- 
strated in Ref.Hniandini 



(31) 



where qi is the net charge in the atomic sphere of 
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the i-th. alloy component, S^s Wigner-Seitz radius, and 
User the screening constant which can be determined 
from the average values of the net charges and elec- 
trostatic potentials of the alloy components in supercell 
calculations.^*' The total energy in this case should 
be also corrected by the energy of the on-site screened 



Coulomb interactions, 
nary alloy is 



which in the case of a bi- 



iS*,, 



(32) 



Here Ci is the concentration of the j-th alloy component, 
and an additional coefficient /3 takes care of the non- 
spherical contributions to the electrostatic energy (/3 = 1 
if the multipole moment contributions to the electrostatic 
energy and potential are neglected)]^ 

To do accurate single-site DFT-CPA calculations of 
the total energy of a random alloy, one should first de- 
termine both screening constants, User and /?. This can 
only be done in supercell calculations, which enables one 
to go beyond the single-site approximation and deter- 
mine the electrostatic energy and potential accurately. 
Obviously, such calculations are computationally very 
demanding, and the screening constants are, therefore, 
usually assumed to be constant for a given alloy, or 
at most have some concentration dependence.'^ In our 
EMTO-CPA calculations here, we determine screening 
constants from the corresponding ELSGF calculations for 
the initial undistorted bcc structure. In the FM calcula- 
tions, we get User = 0.7129 and in the DLM calculations, 
User = 0.778. 

The ELSGF calculations have been done with LIZ=3 
for the undistorted bcc lattice (two nearest-neighbor co- 
ordination shells for every site) and LIZ=5 for the dis- 
torted lattices. The shear moduli have been determined 
from fitting the total energies of the alloy for five distor- 
tions with the step 1 and 0.5 % in the cases of c' and C44, 
respectively. The results of the calculations are presented 
in Table m 



TABLE I. Elastic constants (in Mbar) of the Feo.875Cro.125 
random alloy with and without SRO at the first coordination 
shell obtained by the EMTO-CPA, ELSGF and PAW meth- 
ods in the ferromagnetic and DLM states. 





ELSGF 


EMTO-CPA 


PAW 


Constant 


Qfi = ai — 0.1 


(ai = 0) 


(qi = 0) 


C44-FM 


1.063 1.090 


1.099 


0.900 


C44-DLM 


1.196 1.149 


1.194 




c'-FM 


0.705 0.741 


0.740 


0.555 


c'-DLM 


0.217 0.226 


0.231 





As one can see from the table, the EMTO-CPA cal- 
culations overestimate (in this particular case, of course) 
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FIG. 9. The screening parameter, Q^cr, as a function of the 
distortion parameter in the calculation of c' in the ferromag- 
netic Feo.875Cro.125 random alloy. 



both elastic constants. In the case of ferromagnetic cal- 
culations, this results from the assumption of the inde- 
pendence of the screening constants, ascr on the amount 
of deformation. In Fig. |9) we demonstrate that the as- 
sumption is not really accurate. The screening constants 
have been determined in the corresponding ELSGF cal- 
culations as 



m - V 



(33) 



where (Vi) and (qi) are the average values of the electro- 
static potential and net charge of the atomic sphere of the 
i-th alloy component in the supercell, and V = J2i Ci{Vi) 
{V = in the absence of the multipole-moment contri- 
butions). The screening constant does not depend on 
the alloy component in the case of binary alloys, but 
it becomes component-dependent if the number of alloy 
components is greater than two. 

The screening constant growing with the deformation 
parameter implies the decrease of the screened-Coulomb- 
interaction energy and as a result, of the total energy of 
the alloy. This produces the effect of softening of the elas- 
tic constants. In CPA calculations, this effect is neglected 
and higher values of the shear moduli are obtained. It 
is interesting to note that the phase-separation type of 
the atomic SRO leads to the increase of both shear mod- 
uli by about the same 5% as the error coming from the 
incorrect electrostatics in the EMTO-CPA calculations. 

There is also a noticeable effect of the magnetic state 
on elastic properties. Going over from FM to the DLM 
state in Feo.875Cro.125 results in a significant drop of c' by 
about a factor of 3, and C44-constant increases by about 
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10% (more thoroughly this system is discussed in Ref. 
|41]) . It is clear that the effect of the atomic SRO on the 
elastic constants in FeCr alloys in the paramagnetic state 
is much less pronounced than in the ferromagnetic case. 
The difference between EMTO-CPA and ELSGF calcu- 
lations is also small in the paramagnetic state. This weak 
sensitivity of the strain with respect to the magnetic mo- 
ment, as well as the screening constant, is due to isotropy 
of the system in this case. 

In Table [l) we also show the results of the projector- 
augmented-wave (PAW)^^ calculations of the random 
Feo.875Cro.125 alloy that has been modeled by a 128-atom 
supercell.l^ Clearly, the present implementation of the 
total-energy full-charge-density technique produces quite 
a substantial error, of about 25-30%, and this issue needs 
a separate investigation. Nevertheless, ELSGF seems to 
be a versatile and sufficiently accurate tool to investi- 
gate complex effects of the atomic SRO, especially when 
it comes to the impact of a magnetic state on various 
properties of alloys. 



D. Surface segregation energies of Cr on Fe(OOl) 
and Fe(Oll) 

Apart from direct applications of LSGF to random al- 
loys, the method can also be used to study inhomoge- 
neous systems, such as surfaces, interfaces, impurities 
etc. One of such prospective applications is the calcu- 
lation of solution and segregation energies in cases when 
the size mismatch of alloying components is small (and 
the relaxation energy is small). The advantage of the 
LSGF formalism here is that the LIZ effectively cuts off a 
spurious interaction between alloying species. This prop- 
erty is very important for systems either with strong and 
long-range effective interactions and/or exhibiting strong 
concentration dependence of the alloying behavior. This 
is exactly the case of Fe-Cr alloys, where (as has been 
demonstrated in Ref. |31]), the size of the supercell plays 
a very important role in the determination of the surface 
segregation energy. 

An alternative way to calculate this energy is to use 
the CPA-based method;^ However, in the presence of 
a substantial charge transfer between atomic spheres of 
the alloy components, it is very difficult to accurately 
take care of the on-site screened Coulomb interactions 
for surface alloys, and this can lead to a large error as 
was demonstrated in Ref. 



TABLE II. Surface segregation energies (in eV) of Cr on the 
(100) and (110) surfaces of Fe in the FM and DLM states. 





Fc(llO) 


Fc(lOO) 


FM 


0.068 


0.204 


DLM 


0.144 


0.190 



In this section, we apply the ELSGF method to the 
calculation of the surface segregation energy of Cr onto 
the (001) and (110) surfaces of bcc Fe. Although this 
energy has nothing to do with the corrosion resistance 
of steels, as has been frequently claimed, it is anyway an 
important thermodynamic quantity. 

There exists several calculations of the surface segrega- 
tion energies for this system,^"* 46 49 y^^^ q£ them have 
been done for the ferromagnetic state. An additional ad- 
vantageous feature of LSGF is that the surface segrega- 
tion energy (or more general surface-related quantities) 
can be obtained in the paramagnetic (DLM) state, which 
is important because phase transformations usually take 
place in this state. 

The surface segregation energy is the energy difference 
between two configurations of an impurity atom: one 
with the impurity being in the surface layer and another 
one, when it is in the bulk. The supercell approach is 
then reduced to the calculations of such two systems. The 
surface in this case can be modeled using a slab geometry 
with a vacuum region which, in the case of the EMTO 
method, is filled with empty spheres. In our case, we have 
chosen a 20-layer slab for the (001) surface (13 atomic 
and 7 vacuum layers) and a 14-layer slab (9 atomic and 5 
vacuum layers) for the (110) surface. The corresponding 
20- and 14-atom unit cells have been used to define the 
effective medium, while the entire supercells have been 
constructed from the initial slab unit cells by translations 
in the plane parallel to the surface layer repeating the 
unit cells 6x6 and 8x6 times for the (100) and (110) 
surfaces, respectively. The supercell has consisted, thus, 
of 672 sites for the (110) surface and of 720 sites for the 
(001) surface. 

In the ELSGF calculations, the LIZ has consisted of 
a central atom and its two nearest-neighbor (bcc) coor- 
dination shells (LIZ=3). The self-consistent calculations 
have been done for the room-temperature experimental 
lattice constant of Fe, 2.86 A, using local-density approx- 
imation (LDA).Eal 

The results of the calculations are presented in Ta- 
ble|ll] They can be compared to the results for the (100) 
surface obtained by Ponomareva et al.^ who also did 
calculations for the room-temperature lattice constant of 
Fe using the PAW method and found that the surface seg- 
regation energy of Cr on the (100) surface of Fe is 0.190 
eV. Our ELSGF result, 0.204 eV, is in a good agreement 
with the PAW result. 

The segregation energies obtained in Ref. HHIare signif- 
icantly lower than those in Table [ll) -0.001 and 0.076 eV 
for the (110) and (100) surfaces, respectively, most prob- 
ably because they were obtained for a relatively small 
supercell and for the theoretical equilibrium lattice con- 
stant of Fe. As was demonstrated by Ponomareva et al.^ 
both these parameters strongly affect the results, and we 
therefore believe that our results are quantitatively accu- 
rate. 

Finally, one can see that there is quite a pronounced 
effect of the magnetic state on the segregation energy of 
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the (110) surface: the surface segregation energy in the 
DLM (paramagnetic) state is almost doubled compared 
to that in the FM state. This means that at elevated 
temperatures, relevant to experimentally achievable equi- 
librium, the surface segregation of Cr atoms towards the 
(110) should be reduced. Let us note that this surface 
has an important role in the thermodynamics since it is 
the most closely packed surface in the case of the bcc 
structure, and thus has the lowest energy. It is also clear 
that such a result is extremely difficult, if possible, to ob- 
tain at the same level of accuracy by any other existing 
method. 



IV. CONCLUSION 

We have demonstrated that the ELSGF method can 
be a rather accurate and versatile tool for studying local- 
environment effects in random alloys as well as in inho- 
mogeneous systems, such as surfaces and interfaces. In 
particular, we have applied it to the Fe-rich FeCr alloy 
and found that an intimate coupling between SRO and 
the equilibrium lattice constant as well as elastic prop- 
erties results from the strong sensitivity of Cr magnetic 
moments on the local environment. Among the advan- 
tages of the method is its order- scaling, which makes 
the implementation easy to parallelize and enables one to 
treat large supercells consisting of A^ '--^ 10* — 10^ atoms. 
In addition, the capability to treat the high-temperature 
paramagnetic state renders possible investigating phase 
transitions in magnetic alloys. 

Compared to the KKR-ASA implementation, ELSGF 
has a much more accurate normalization of states and 



can, therefore, be applied to systems with a distorted 
structure and large ion-size mismatches. This feature 
opens up a completely new possibility to take into ac- 
count random relaxations in alloys within ELSGF. For- 
mally, this amounts to introducing an additional pertur- 
bation of the structure constants into the Dyson equa- 
tion, Eq. (18). Accurate evaluation of the kinetic en- 



ergy is a necessary prerequisite for such an expansion 
in the structure-constant perturbation to give good re- 
sults, which can be achieved by a full-charge-density-self- 
consistent implementation of ELSGF. 

Other implementations of LSGF are possible. For 
instance, a fully-relativistic EMTO-LSGF method can 
be realized by replacing Green's functions in spin-up- 
spin-down channels with full Dirac 2 x 2-spinor Green's 
functions. A fuU-potential-KKR implementation is also 
straightforward. 
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